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Abstract 

A recently proposed regularization-independent method is used for the first 
time to solve the renormalized fermion Schwinger-Dyson equation numeri- 
cally in quenched QED4. The Curtis-Pennington vertex is used to illustrate 
the technique and to facilitate comparison with previous calculations which 
used the alternative regularization schemes of modified ultraviolet cut-off and 
dimensional regularization. Our new results are in excellent numerical agree- 
ment with these, and so we can now conclude with confidence that there is no 
residual regularization dependence in these results. Moreover, from a compu- 
tational point of view the regularization independent method has enormous 
advantages, since all integrals are absolutely convergent by construction, and 
so do not mix small and arbitrarily large momentum scales. We analytically 
predict power law behaviour in the asymptotic region, which is confirmed nu- 
merically with high precision. The successful demonstration of this efficient 
new technique opens the way for studies of unquenched QED to be undertaken 
in the near future. 
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I. INTRODUCTION 



The divergences inherent in quantum field theories have plagued physicists for years. 
The infinities are removed in a two-step process; first the divergences are controlled by 
a regulator, then the regulator is removed using the renormalization procedure to obtain 
finite, renormalization-independent physical quantities. Depending on the type of regulator 
introduced, different problems can occur. It is useful to discuss the difficulties that arise 
in the Schwinger-Dyson equations (SDE's) since these are identical to the difficulties 
encountered in perturbation theory and because we will be using them as a tool to study 
non-perturbative QED. In the literature SDE's in quenched QED 4 are commonly studied 
by using the ultraviolet cut-off regularization scheme 0. One difficulty faced is that the 
use of an ultra-violet (UV) cut-off A to regulate the integration will in general lead to 
an explicit violation of gauge covariance 0. Because this regularization scheme does not 
respect translation invariance in the loop-momentum integration, it will lead to an explicit 
gauge-covariant violating contribution in the result, even after A is taken to infinity. More 
precisely, this violation of gauge invariance has been observed in quenched QED calculations 
employing the Curtis-Pennington (CP) |14| electron-photon vertex and may be traced back 
to a certain, logarithmically divergent, 4-dimensional momentum integral which vanishes 
because of rotational symmetry at all A < oo, but leads to a finite contribution for A — > oo. 
It is this discontinuous behaviour as a function of A which complicates correct numerical 
renormalization with this regulator. Incorrect results will be obtained unless care is taken 
to identify and remove 'gauge-covariance violating terms' 0- [[F]. In its favour, cut-off 
regularization is computationally economical compared to dimensional regularization, and 
gives accurate answers after gauge-covariant violating terms are removed. 

On the other hand, SDE studies implemented using a gauge-invariant regularization 
scheme, such as dimensional regularization ||, will not have such a problem. However it 
will cause the generation of dynamical mass for all coupling constants in D ^ 4, instead 
of only above critical coupling in D = 4. In non-perturbative studies this scheme is com- 
putationally more demanding since a careful and time-consuming removal of the regulator 
must be performed involving an extrapolation of many high-precision solutions for different 
e to e = 0. Additionally the accuracy of the results is, in practice, more limited because the 
integration range neccessarily extends to infinity. 

In this paper we will be employing for the first time the regulator-independent method 
recently proposed by Kizilersii et.al. ||. This method deals with the renormalized quantities 
only, as the regulator is removed analytically. The dependence on the mass scale introduced 
by the regulator is traded for the momentum scale /i at which the theory is renormalized 
before performing any numerical calculations. In this way dimensional regularization (or 
any other regulator which doesn't violate gauge covariance) can be used and the regulator 
removed analytically before any numerical calculations are begun. More importantly from a 
numerical pointof view, removal of the regulator means that no longer does one have to solve 
integral equations involving mass scales of vastly different orders of magnitude (and then, 
in addition, take a limit in which one of the scales goes to infinity). Rather, the important 
scales in the problem become scales of physical importance, such as the renormalized mass 
m M and the renormalization scale fi at which this mass is defined. It is therefore to be 
expected that the dominant contributions to any integrands will be from a finite region of 
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momenta. This feature is generic; i.e., it is independent of the particular vertex ansatz that 
one makes use of, and it remains a valid consideration for an arbitrary renormalizable field 
theory. 

The regularization-independent approach is computationally very economical and very 
accurate; however the price paid is that one loses all contact with the bare theory. In particu- 
lar, because of this, one cannot study dynamical chiral symmetry breaking in this approach 
by simply setting the bare mass mo = and investigating at what value of coupling the 
dynamical fermion mass is generated. It is in this crucial point we differ from the analysis 
in Ref. |TD[ where, within quenched QED, a removal of the regulator was attempted in a 



manner which has some similarity to what we do here. Indeed, as emphasized by Miransky 
(Refs. fill , jl2| as well as chapter 10.7 of Ref. [13|]), some care needs to be taken with the 



treatment of the bare mass while removing the regulator in order to avoid drawing incorrect 
conclusions about the presence or absence of dynamical chiral symmetry breaking in gauge 
theories. In the present approach this problem is avoided. Dynamical chiral symmetry 
breaking, both explicit as well as dynamical, is characterized by a non-zero renormalized 
mass m M , and hence this quantity in itself cannot distinguish between these two possibili- 
ties. In section |III A| , we use (within quenched QED4) the appearance of oscillations in the 



renormalized fermion mass function [|TT| as the indicator of the onset of dynamical symme- 
try breaking. All previous studies have confirmed that the occurrence of dynamical mass 
generation and UV oscillations coincide in quenched QED 4 . 

In this paper we demonstrate this approach numerically for the specific case of quenched 
QED4 using the CP vertex to test its validity against previous works. Our goal is to ulti- 
mately implement this technique in studies of unquenched QED4, which have now become 
numerically tractable because of this approach. 

In section |TI] we formulate the regularization-independent method for renormalized SDE 
for quenched QED in 4-dimensions with an arbitrary covariant gauge. Section |T| discusses 
the large momentum behaviour of the fermion propagator and gives its analytical form. 
We go on to solve SDE's numerically in Euclidean space for the fermion wave-function 
renormalization function and the mass function employing the CP vertex. The integration 
range in these equations is taken to infinity analytically and subsequently evaluated using 
a combination of numerical and analytic results. We do this extrapolation to infinity to 
facilitate comparison with previous results at very high momentum scales and to compare 
these results with their asymptotic forms. We demonstrate numerically to high precision 
the agreement between the regularization independent approach and the modified cut-off 
and dimensional regularization schemes. Finally, in section [TV] we summarize our results 
and conclude. 

II. REGULARIZATION-FREE FORMALISM IN QUENCHED QED 4 

The renormalized Schwinger- Dyson equation for the electron propagator (Fig. [I]) can be 
formulated as 



ZMS°-\p) - iZ 1 ( f i)E(p) . (2.1) 
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FIG. 1. The Schwinger-Dyson equation for the fermion propagator. 



Here Z\{\x) and Z 2 (fi) are the vertex and fermion wave-function renormalization constants 
respectively. These renormalization constants relate the regularized but unrenormalized (i.e. 
bare) and renormalized propagator and vertex by 

r„(*,p;/*) = ZiOi^rJ-topjA) , (2.2) 
S(p;/x) = ^A)S ta (p;A) . (2.3) 

In quenched QED, there is no renormalization of the electron charge (e§ = e 2 = 47ra) and 
the appropriate photon propagator is simply the tree-level form 



DIM = ~J2 



9a0 9— + ? n - 



(2.4) 



where £ is the covariant gauge parameter. We use S(fi; p) to denote the full fermion prop- 
agator renormalized at the momentum scale fi. It can be expressed in terms of two scalar 
functions, Z([i 2 ;p 2 ), the fermion wave-function renormalization function, and M(p 2 ), the 
mass function, by 

s( ^-0rm ■ < 2 - 5 > 

Note that S°(p) = — mo) is the tree-level fermion propagator (i.e., the bare fermion 
propagator in the absence of interactions). The full (proper) renormalized fermion-photon 
vertex is r M (/x; k,p). Multiplying the fermion SDE Eq. ( |2.1| ) by $ and / and taking their 
respective spinor traces, we can separate the fermion self-energy X(p 2 ) into Dirac odd, S^(p 2 ), 
and Dirac even, £ s (p 2 ), parts; i.e. S(p) =p l Sd(p 2 ) + S s (p 2 ). This gives 

V; p 2 ) = Hp) - Hp) W) , (2.6) 

M(p 2 ) Z-\^-p 2 ) = Z 2 {n) m + Z x {p) S s (p 2 ) , (2.7) 

and these self energies make use of any regularization scheme which does not violate gauge 
covariance. Here, the bar over these quantities indicates that we have explicitly separated 
out the renormalization constants Zi(fi); note that for notational brevity we do not indicate 
the implicit dependence on /i 2 of S ds (p 2 ) through the function Z(fi 2 ;p 2 ). 

Because of the Ward-Takahashi identity, which must be satisfied by any acceptable vertex 
Ansatz for T A1 (/i; k,p) and any acceptable renormalization scheme, one has Zi(fi) = Z 2 (fi). 
Making use of this fact, Eqs. ( |2.6| ) and ( |2.7| ) can be rearranged as 
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M(p 2 ) = m + 



)-Z(/i 2 ;p 2 )S (i (p 2 ) , 
M(p 2 )E d (p 2 ) + £ s (p 2 ) 



(2.8) 
(2.9) 



In order to avoid cumbersome notation, we have not explicitly indicated functional depen- 
dence on the regulator in Eqs. (|2.1|) - (|2l3D but it should be understood. The renormalization 
constants Z 1)2 (/x) and the bare fermion mass m§ are regulator dependent in the above equa- 
tions. As one removes the regulator, the integrals on the right hand side of Eq. (|2.1|) , and 
hence H d (p 2 ) and Y> s (p 2 ), diverge logarithmically. It is the defining feature of a renormal- 
izable field theory that these divergences may be absorbed into the constants Z 12 (/i) and 
into the bare mass mo, rendering finite, regularization-independent limits for Z(/j, 2 ;p 2 ) and 
M(p 2 ). However, we can make use of the p 2 independence of Zi^ifi) and mo in order to 
eliminate these constants from the above equations. The renormalization conditions for the 
fermion wave-function renormalization function and the mass function are 



Z(/i 2 ;/i 2 



1 and M(/i 2 ) = m p 



Evaluating Eqs. (|2.8|) and (|2.9| ) at a second momentum which, we take to be p 2 
taking the difference one obtains the central equations that we solve here 



Z^ 2 -p 2 ) 


= 1 + Z(p 2 ;p 2 )Z d (p 2 ) - S d (/i 2 ) 
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M(p 2 ) 


= m M + 


'M{p 2 )Y, d {p 2 ) + E s (p 2 ) 




m M S d (/i 2 ) + £ s (/i 2 ) 





(2.10) 
/i 2 , and 



(2.11) 



As the left hand sides of these equations must be finite as the regulator is removed, then the 
right hand side must be also, even though the individual terms on the RHS may separately 
diverge. 

These renormalized equations Eq. (|2.11| ) with the regulator removed provide a starting 
point for non-perturbative investigations which have significant advantages over the usual 
treatment found in the literature @- ||. In the following section we shall illustrate this 
approach by turning to the example of quenched QED with the Curtis-Pennington vertex. 



III. QUENCHED QED WITH THE CURTIS-PENNINGTON VERTEX 

In this section we use the regulator-independent method for solving the renormalized SDE 
in Eq. ( |2.11| ) in an arbitrary covariant gauge. As usual we write the full vertex as a sum of 
the Ball-Chiu vertex [15[ (longitudinal part satisfying the Ward-Takahashi identity) and the 
Curtis-Pennington term |14j] (the transverse part satisfying multiplicative renormalizability 
of the fermion propagator), 

r»(wk,p) = ri iC (wk,p) + Te(wk,p) [^( P 2 -k 2 ) + (p + kr^- ff)] , (3.i) 

where r^ c .(/i; k,p) and T 6 (fi;k,p) are given in appendix [A]. After Wick rotating into Eu- 
clidean space and performing the angular integration, Eq. ( p.ll| ) can be written as 
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i - * r de 1 

4vr J P 2 [k 2 + M 2 {k 2 )} 



a foe 
+ — 



dk 2 



a 



M(p 2 ) = m M + — 



An Jo [k 2 + M 2 (k 2 )} 
dk 2 



An Jo [k 2 + M 2 (k 2 



Z(j,V)/(*V) - /(*V)] , (3.2) 

J(A; 2 ,p 2 )- J(A: 2 ,/i 2 ) 
+ M(p 2 )I(k 2 ,p 2 )- mil I(k 2 ,^)} , (3.3) 



where the kernel functions I(p 2 , k 2 ) and J(p 2 , k 2 ) are also given in appendix [A]. The kernels 
have the same form as Ref. |16j except that the gauge-covariance violating term has been 
removed, since it does not survive the four-dimensional momentum integration in the absence 
of a cut-off. This term does not vanish if the angular integral is done first with the radial 
integral taken to infinity afterwards, which is why it must be removed by hand in the cut-off 
approach. 

The above equations are finite which can be seen by analysis of the large momentum 
limits of the integrands for Z and M, which behave asymptotically as (k 2 )~ n and (k 2 )~ r 
respectively, where n,r > 1. The equation for the wave-function renormalization function 
Eq. ( |3.2| ) is the same as Eq. (10) in Ref. [ICfl, however our treatment of the mass function 
differs from theirs. 



A. Asymptotic limits of the solutions 

In this subsection we solve Eqs. ( |3.2p and ( |3.3| ) analytically for momenta p 2 much larger 
than M 2 (p 2 ) (the asymptotic region) by linearizing them for finite mass, in a manner similar 
to that used by Atkinson et al. [|l(| for the chirally symmetric theory (i.e. m = 0). It is 
convenient to temporarily take the renormalization point \x to also be very large, removing 
this constraint once our derivation is complete. The integrals in Eqs. ( |3.2D and (|3.3|) are 
dominated by contributions around k 2 = p 2 and k 2 = /i 2 , so when p 2 and fi 2 are in the 
asymptotic region k 2 is neccessarily also much greater than M 2 (k 2 ). Expanding in powers 
of M 2 {k 2 )/k 2 and keeping at most linear terms in M(k 2 ), the kernel functions take their 
asymptotic forms I(p 2 ,k 2 ) — > and J(p 2 ,k 2 ) — > J'(p 2 ,k 2 ), where J'(k 2 ;p 2 ) is defined in 
appendix |A[ Hence when p 2 and /i 2 are both in the asymptotic region, we obtain 

z^ 2 ;p 2 ) = i - ^z^k 2 ) , (3.4) 

M{p 2 ) = m M + — / dk 2 \j\k 2 ,p 2 ) - J\k 2 ,fi 2 )} . (3.5) 

47T JO L J 

These linearized equations are scale invariant and admit power law solutions. 
The solution of this asymptotic form of the Z equation is easily seen to be 

Z^ 2 -p 2 ) = [4) with v=^ , (3.6) 



which differs from Ref. [16] due to the absence of the gauge covariance violating term. We 



thus obtain for p 2 and p both large that 
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w> (ELY with „ = £ , (3 ,) 



Z(/i 2 ;p' 2 ) Vp'V 47r 

which is valid for an arbitrary renormalization point fi since the ratio is by definition 
renormalization-point independent. It then follows that for large p 2 and any arbitrary /i 2 
that 

Z^ 2 ;p 2 ) = with u=^- , (3.8) 

for some appropriate such that — > 1 as /i 2 enters the asymptotic region. 

We shall see that power law solutions of the M equation occur in two regimes, depending 
on the value of a. Those with a real exponent can be identified with the subcritical regime, 
where a is less than its critical value a c which marks the onset of dynamical chiral symmetry 
breaking. Those with complex exponent correspond to oscillating solutions which correspond 
to the supercritical regime where a > a c . We examine these two cases separately in what 
follows. 



1. Subcritical asymptotic solution 



Where, as before, /z 2 and p 2 are understood to be chosen much greater than all other 
mass scales, we try a solution of Eq. (|3.5| ) of the form 



M{p 2 ) 



m. 



p 



seR. 



(3.9) 



Solving for s, we find 

27T COt S7T — 7T COt h>7T + 37T COtfl^ — sW + 

2£ 



(1-s) + 



+ 



V 



+ 



s-1 
u-s + 1 







(3.10) 



This equation is the same as in Ref. but in contrast to the unrenormalized equation 
examined in that paper, this power Ansatz is valid for the renormalized equation that we 
have here whether or not the bare mass m is zero. We can again form a ratio, i.e., provided 
p 2 and p' 2 are much larger than all other mass scales then we have 



M(p 2 ) 
M(p 12 ) 



p" 



n 



seR 



(3.11) 



Both Eqs. ( p.lO| ) and ( 3.11|) are manifestly renormalization-point independent and so are 
valid for an arbitrary choice of renormalization point \x. It then follows that for large p 2 and 
any arbitrary /i 2 that 



M(p 2 



D, 



p 

J 2 



s,D^e R 



(3.12) 
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for some appropriate D^. We see that D^Kn 2 ) -8 is a renormalization-point independent 
constant and that we must have — *> m M as /i 2 moves into the asymptotic region. 

There is in fact more than one real solution to s in Eq. ( |3.10| ) (see also Ref. ||16||), e.g., 



in Landau gauge there are two. However, only one of these matches smoothly onto the 
perturbative solution, and hence is the relevant solution. The other solutions appear as a 
spurious byproduct of the linearized approximation in Eq. ( |3.5|) (Solutions which do not 
match smoothly onto perturbation theory can arise in Eq. ( ^.5[ ) if the integrals diverge like 
l/a, due to infrared divergences; this cannot occur in Eq. ( |3.3| ) because of the regulating 
mass in the denominator). Since s is real there are no oscillations in the mass function and 
hence there is no chiral symmetry breaking for these couplings. While Atkinson et. al. are 
unable to generate a mass in this region, due to their initial condition being tuq = 0, we are 
able to do so since we are not confined by this constraint. 



2. Supercritical asymptotic solution 

In the case where s is complex, the (real) mass function is a superposition of powers of 
s and its complex conjugate. For large p 2 and arbitrary jj? we have as the counterpart to 



Eq. (£TJ) 




M(p 2 ) = -D, ( ^ ) • s,D,EC (3.13) 

where is some appropriate complex constant, and 

l -(D, + D;) = Re(£>„)-ro„ (3.14) 

as /i 2 enters the asymptotic region. Here s and s* are complex conjugate solutions to 
Eq. ( |3.10| ) and the combination results in a mass function that oscillates. This oscillatory 



behaviour is more transparent if we write Eq. ( |3.13| ) as 

/ 2 \ -Re(s) 

M (P 2 ) = ( J { Re ( D J cos [ Im ( s ) lo S (V/V)] + Im P M ) sin [im (s) log (p 2 /fi 2 )} } 

(3.15) 

which has the form of a phase-shifted cosine function periodic in log(p 2 /// 2 ), with a period 
of 2 7r/Im(s), modulated by an exponentially decaying envelope of (p 2 /yU 2 ) -Re ^. The oscil- 
lations are an indication that chiral symmetry is broken in quenched QED employing the 
CP vertex at sufficiently strong coupling a > a c . 



B. Numerical solutions 

In the regularization-independent approach, the equations are constructed so as to re- 
move all infinities from the outset and so the UV region of the integral is not a crucial 
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limiting factor. Moreover, establishing exact asymptotic forms of the wave-function renor- 
malization and mass functions in the previous subsection make it possible to analytically 
integrate Eqs. ( |3.2| ) and ( |3.3| ) from the highest grid point to infinity. We denote this highest 
grid point, where the asymptotic analytic forms are matched onto the numerical results, as 
the "matching point" k^. Hence we can rewrite Eqs. ( j3.2|) and ( |3.3| ) as 



where the matching point, k^, is chosen sufficiently large that the asymptotic formulas are 
valid (k 2 m ^> M 2 (/c^)). These equations serve as definitions of the analytic forms that we need 
to evaluate, i.e., Z high (/z 2 ; k^p 2 ) and M high (/^ 2 ; k^p 2 ) are defined as the contributions arising 
from integrating from k 2 m to infinity for Eqs. ( |3.16| ) and ( |3.17| ) respectively. The analytic 
forms for Zhigh an d -^high an d their derivations are given in appendix |B]. By calculating 
Z(n 2 ;p 2 ) and M(p 2 ) in this way we can achieve very high accuracy in the UV region. 

Eqs. ( [3.16D and (|3.17|) have been solved numerically in Euclidean space for Z(fj, 2 ;p 2 ) and 
M(p 2 ) with a variety of gauges £, renormalization points fi 2 and renormalized masses m M for 
the couplings a = 0.6 (subcritical) and a = 1.5 (supercritical). Each solution was iterated 
from an initial guess until Z(fi 2 ;p 2 ) and M(p 2 ) converge; our convergence criteria is that the 
change in successive solutions is less than one part in 10 8 at each momentum point. Every 
iteration the data is refitted to the asymptotic analytic forms of Z(fi 2 ;p 2 ) and M(p 2 ). 

We have found, in line with the above discussion, that the regularization-independent 
method allows us to extend solutions for momenta up to 10 65 as opposed to solutions obtained 
from cut-off regularization where numerical round-off error meant that the momenta in our 
solutions could typically not exceed O(10 18 ) 0, |>|, 0. We also verified that our solutions 
do not change when we vary the location of the point k m where the matching of our numerical 
and analytical solutions is carried out. 

In Figs. and |3|, we show typical solutions, using regularization independent regular- 
ization, of the fermion wave-function renormalization and mass functions for subcritical 
(a = 0.6) and supercritical (a = 1.5) cases respectively. Also shown are corresponding 
solutions based on the theoretical asymptotic powers with fitted scales. The log-log nature 
of the figures emphasizes the UV behaviour of the propagator. One can see there is excellent 
agreement in the region p 2 3> M 2 (p 2 ). To emphasize this point we present the theoretical 
and numerically calculated powers for these two solutions in Table |. 

In summary : our knowledge of the exact forms of Z(/i 2 ;p 2 ) and M(p 2 ) in the asymptotic 
region saves us from the need to extend our numerical solutions far into the ultra-violet. 
Therefore Eqs. ( |3.16| ) and ( |3.17| ) are solved numerically up to the matching point {k 2 ^ 
in the UV, whereafter the analytical solution joins smoothly to the numerical one. The 




+ Z Klgh {n 2 ]k 2 ml p 2 ), 



(3.16) 




+ M{p 2 ) I(k 2 ,p 2 ) - m, I(k 2 ,fi 2 ) + M high Gu 2 ; k 2 m ,p 2 ) 



(3.17) 
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fitting parameters [C^, Re (Dp), Im (D^), Re (s), Im (s), u] for the analytic continuation are 
re-calculated after each iteration. 
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FIG. 2. Comparison of the numerical solution for the finite renormalization 1 / Z \ji 2 ; p 2 ) , and 
the mass function M(p 2 ) (solid lines) and their predicted asymptotic behaviour with matching 
scales (dashed lines) from Eqs. ( |3.8| ) and fl3.12|) for the subcritical coupling a = 0.6. The example 
solution was for a renormalized mass = 400 (arbitrary units), renormalization point /i 2 = 10 8 
and gauge parameter £ = 0.25. 
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(dashed lines) from Eqs. ( |3,8| ) and ( 3.13j ) for the supercritical coupling a = 1.5. The example 



solution was for a renormalized mass m« = 400 (arbitrary units), renormalization point fi 2 = 10 8 
and gauge parameter £ = 0.25. The " guide lines" are added simply to guide the eye between data 
points on this logarithmic plot 
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TABLE I. This table compares the powers of the asymptotes of the finite renormalization and 
mass functions in Figs. || and |3[ as determined analytically and by fitting the numerical solutions 
to a power law form. 



Figure 


coupling 


determination 


V 


Re(s) 


Im (s) 


1 


a = 0.6 


analytical 
numerical 


0.01193662073 
0.01193662073 


-0.181667015 
-0.181666808 




I 


a = 1.5 


analytical 
numerical 


0.02984155182 
0.02984155183 


-0.416012578 
-0.416012521 


0.418128942 
0.418129001 



C. Comparison of regularization schemes 

We can now compare numerical solutions from the regularization-independent approach 
to those from cut-off regularization with the gauge-covariance modification @, 0, @ and 
those from dimensional regularization M. In cut-off regularization the fermion self-energies 
are integrated on a logarithmically spaced grid in k 2 momentum up to the highest (cut- 
off) momentum A 2 . Additionally in dimensional regularization, an estimate is made of the 
contribution to the integral from the highest grid point to infinity. In both the modified 
cut-off and the dimensional regularization studies subtractive renormalization is performed 
numerically for the regularized (but otherwise divergent) fermion self-energies. For the 
modified UV cut-off regularization approach it was necessary to perform the calculation 
at several values of A and in principle perform a A — > oo extrapolation. In practice it 
was found that it was sufficient simply to ensure that A was chosen large enough. For the 
dimensional regularization approach it was necessary to calculate solutions at high accuracy 
for very many values of e and then carefully extrapolate e to 0. Numerical limitations made 
it difficult to obtain solutions at very small e, which in turn limit the achievable accuracy of 
the e — > extrapolation. 

In Figs. |] and [5], we compare these three regularization methods for subcritical (a = 0.6) 
and supercritical (a = 1.5) cases respectively, with the standard parameter choice of £ = 
0.25, /i 2 = 10 8 , m M = 400. One can see the agreement between the regulator-independent 
(NR) and the modified cut-off and dimensionally regularized solutions is excellent. They 
are indistinguishable on the main figures: the inserts in Figs. f| and [5] have the same p 2 scale 
and reveal the remarkable agreement between them in the infrared region. In Tables [TT] and 



ITT| we quantify the relative differences achieved between the three regularization schemes at 
a variety of momentum values. The difference between the results is entirely attributable to 
the limitations achievable in numerical precision. We can now conclude with some confidence 
that the three regularization schemes give identical results for the renormalized solutions. 
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Subcritical Regularization Comparison 
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FIG. 4. The finite renormalization l/Z(fi 2 ;p 2 ) and dynamical mass M(p 2 ) for the solution of 
the fermion SDE for subcritical coupling a = 0.6 and gauge parameter £ = 0.25 found from the 
regularization-independent (NR) method compared with solutions using the modified UV cut-off 
regulator and dimensional regularization. The dimensional regularization solution shown is the 
result of extrapolating various finite e solutions at scale 10 3 to e = using a fit cubic in e at 
each momentum point. All solutions have renormalized mass = 400 (in arbitrary units) at the 
renormalization point fj? = 10 8 . The small variation between the three regularization schemes is 
entirely attributable to limitations achievable in numerical precision. 
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Supercritical Regularization Comparison 
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FIG. 5. The finite renormalization l/Z(fi 2 ;p 2 ) and dynamical mass M(p 2 ) for the solution of 
the fermion SDE for supercritical coupling a = 1.5 and gauge parameter £ = 0.25 found from the 
regularization-independent (NR) method compared with solutions using the modified UV cut-off 
regulator and dimensional regularization. The dimensional regularization solution shown is the 
result of extrapolating various finite e solutions at scale 10 3 to e = using a fit cubic in e at 
each momentum point. All solutions have renormalized mass m M = 400 (in arbitrary units) at the 
renormalization point /i 2 = 10 8 . The small variation between the three regularization schemes is 
entirely attributable to limitations achievable in numerical precision. 
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TABLE II. This table shows the relative differences achieved for the finite renormalization 
1 / Z (/j, 2 ; p 2 ) between the regularization-independent approach and the modified UV cut-off and 
dimensional regularization (for cubic and quartic fits) approaches. The differences are averages 
over all momentum points solved in each order of magnitude shown. The solution parameters are 
the same as those in Figs. ||| and [|. 
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TABLE III. This table shows the relative differences achieved for the mass function M{p 2 ) 
between the regularization-independent approach and the modified UV cut-off and dimensional 
regularization approaches. The parameters are the same as those in Figs. || and |5L 
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IV. CONCLUSIONS AND OUTLOOK 

In this paper we have for the first time solved the Schwinger-Dyson equations for the 
fermion propagator in quenched QED 4 using the regularization-independent approach re- 
cently proposed in Ref. 0. This has been done for the particular choice of the Curtis- 
Pennington transverse photon-fermion vertex, since this facilitates comparison with previous 
results which used the (gauge-covariance) modified UV cut-off and the dimensional regular- 
ization schemes. We have carried out precise calculations in these three approaches and have 
achieved excellent numerical agreement between them. This clearly demonstrates that we 
are able to achieve high-precision non-perturbative calculations of the renormalized fermion 
propagator that are free from any spurious errors which might arise from the regularization 
procedure itself. 

We have derived and used the asymptotic analytic form of the solutions to obtain high ac- 
curacy even at extremely large momentum scales (O(10 65 )). The reason that this is possible 
now is because in the regularization-independent approach all momentum integrations are 
finite by construction, [see Eq. ( [2.11 )1 . No bare mass or renormalization constants appear 
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in this formulation, since they have been eliminated by combining and subracting renormal- 
ized quantities. Since contact with the bare theory is lost, the onset of dynamical chiral 
symmetry breaking is signalled by the onset of oscillations in the UV mass function. This is 
a well-studied phenomena in quenched QED 4 . We have derived the explicit analytical forms 
for the oscillations in the asymptotic region above critical coupling, including the period and 
decay envelope of these. 

The importance of this regularization-independent approach lies in the fact that since 
all unregularized momentum integrations are finite from the outset, we do not have the 
mixing of small and arbitrarily large momentum scales in the intermadiate stages of our 
numerical calculations. This means that we can achieve high accuracy for solutions in the 
low and medium momentum regime with great numerical economy. This new approach 
will now permit numerically tractable studies of unquenched QED4. These studies are now 
underway. 
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APPENDIX A: 



The expressions mentioned in section |TTT] are given below. 
The Ball-Chiu vertex [15 



is 



r^ c (/x;fc,p) 



1 1 

+ 



2 \Z(^i 2 ;k 2 ) Z(/z 2 ;p 2 ) 



r 



(k+py 
k 2 — p 2 



+ f) ( M(k 2 ) M(p 2 ) 



Z(/2 2 ;k 2 ) Z(n 2 ;p 2 )J 2 \Z(fi 2 ;k 2 ) Z(/j, 2 ;p 2 ) 

The coefficient function of the transverse vertex (Curtis-Pennington) Jl4| is 

1 \ 



T 6 (n;k,p) = -— 



2d\Z{p?-k 2 ) Z(n 2 ;p 2 ) 



where 



d 



{(P _ p 2f + ( M 2 (A; 2) + M 2 (p 2 )f) 
(k 2 + P 2 ) 



(Al) 



(A2) 



The constituents of the integrand in Eqs. ( |3.2|) and (|3.3|) are the renormalization-point 
independent kernel functions 
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M{k 2 



+ 



2 (k 2 - p 2 

1 (k 2 + p 2 )(M 2 (k 2 ) + M 2 (p 2 )) 2 
2{(k 2 - p 2 ) 2 + (M 2 (k 2 ) + M 2 (p 2 )) 2 } 



M(k 2 )-M(p 2 ) Z ^ k ^ 



Z{\i 2 \p 2 ) 

Z{^k 2 ) 
Z{jj?;p 2 ) 



9{p 2 - k 2 ) + 6{k 2 - p 1 



Z(fi 2 ;p 2 



k 2 



p^ 



(A3) 



J(k 2 v 2 ) = 3 - M(k 2 ) 1 1 + Z ^ k2) + ( fc4 ~^ 4 ) d _ 

y ,P } 2 1 J \ Z{fi 2 ;p 2 ) {{k 2 -p 2 ) 2 + {M 2 {k 2 )+M 2 {p 2 )) 2 }\ Z(fi 2 ;p 2 ) 



3 2 Z(fi 2 ;k 2 )M(k 2 )- M(p 2 )) k 4 



T Z(/2 2 ;p 2 
Z{^k 2 ) 



M{k l 



k 2 — p 2 
(k 2 



p* 



6{p 2 - k 2 ) + 6{k 2 - p 2 ) 



p z 



0{p 



2 k 2 ) 



Z(fi 2 ;p 2 ) 

Linearizing Eqs. (|A3|) and (|A4fj in terms of the mass function yields 

J (k 2 ,p 2 ) k ^ M , l'(k 2 ,p 2 )=0 , 

J(k 2 ,p 2 ) 

where 

J'(k 2 ,p 2 ) - 



Z i I^Sv M(k 2 ) 6(p 2 - k 2 ) 



+ 



p 2 Z(/j 2 ;p 2 ) 

3(2M(k 2 ) r 2 Z(/2 2 ;k 2 ) 
2 i (p 2 - k 2 ) P Z(fi 2 ;p 2 ) 



- k 



d(p 2 - k 2 ) d(k 2 - p 



+ 



k 2 



Z(fi 2 ; k 2 ) M(p 2 ) - M(k 2 
Z(fi 2 ;p 2 ) (p 2 — k 2 ) 



- 2 0{p 2 -k 2 ) + V ^6{k 2 -p 2 ) 

pZ fcZ 



(A4) 



(A5) 



APPENDIX B: 

This appendix is devoted to the analytic calculation of the wave-function renormalization 
and mass functions in the asymptotic limit. The quantities that we need to evaluate are 

Z high ^ 2 ;k 2 m y) = Q^VMW) [ziM^P 2 )- 1 ^^ 2 )] 
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M high (/i 2 ;/c^,p 2 ) = 



1 



a 



dk 2 



M(Jfe2j47T Jkl k 2 + M 2 (k 2 ) 

J(k 2 lP 2 ) - J(k 2 , /x 2 ) + M(p 2 ) I(k 2 ,p 2 ) - m^I{k 2 , /i 2 ) 



(Bl) 



where I(k 2 ,p 2 ) and J(k 2 ,p 2 ) are given in appendix Provided A; 2 ^ is sufficiently large that 
^> M 2 (/c^) then these quantities become 
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/ 

■7 A' 



k?n k 2 + M 2 (k 2 



j(k 2 , P 2 : 
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(B2) 



We have already shown for k 2 ^> M 2 that the wave-function renormalization and the mass 
function have a power law behaviour 



Z(^k 2 ) = cJ^J 



1 



1 



' k 



2\- s 



m 2 ) = -DA-A +m - 



2~ ^ \ yU 2 



2"^ U 2 



(B3) 



The results for Z high and M high for large fc 2 ^ and arbitrary p 2 and /i 2 can then be given in 
terms of hypergeometric functions 



with 
I(P 2 ) 



Z h ^{n 2 ]k 2 m ,p 2 ) = 
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(B4) 



* J /;,2 \ I 
m/ [ ' 



(B5) 
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D, 



F(l,l + s,2 + s,p 2 /k 2 
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